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A summary of new experimental results and recent theoretical developments discussed 
in the "Hadronic Final States and QCD" working group is presented. 

1 Introduction 

The study of hadronic final states at colliders provides a broad spectrum of insights on 
hadronic structure and spectroscopy and on the physics of quantum chromodynamics (QCD). 
Especially for lcpton-hadron and hadron-hadron colliders, these results are crucial for the 
successful interpretation of any kind of measurement, since QCD effects are omnipresent 
due to the hadronic initial state. In this working group, a variety of new results on many 
different aspects of hadronic final states were presented. These arc summarised in this talk. 



2 Vector Boson Production 

The measurement and modelling of inclu- 
sive single vector boson + jets production 
is both a test of pcrturbativc QCD and 
is an important background measurement 
for new physics. Tevatron measurements 
of Z/-f* + jets and W + jets were pre- 
sented [T|. Figure Q] compares the measured 
inclusive jet cross section from events with 
a Z/'j* — ► e + e~ to the next-to-leading or- 
der (NLO) QCD prediction (open circles). 
The measurement for Nj e t > 1 is scaled up 
(factor 20) for clarity. The shaded bands 
show the total systematic uncertainty ex- 
cept for the 5.8% luminosity uncertainty. 
The dashed and dotted lines indicate the 
uncertainty due to parton distribution func- 
tions (PDFs) and the scale variations of the 
NLO QCD predictions respectively. 

The LHC will be a W/Z factory, pro- 
ducing roughly 20 million VF's and 2 million 
Z's visible by the general purpose detectors 
in each fb _1 of data. The ATLAS exper- 
iment at the LHC has done detailed stud- 
ies [2] using both PYTHIA and ALPGEN 
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Figure 1: (a) Measured inclusive jet cross sec- 
tion as a function of pt in Z /j* — > e + e~ + jets 
events. Plots (b) and (c) show the data/theory 
ratio. Figure taken from pQ. 
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Monte Carlo programs to extract predictions on reconstruction and trigger efficiency, back- 
ground contributions and systematic uncertainties for W/Z + jets events. The dominant 
experimental uncertainty comes from the jet energy scale, and the ultimate ATLAS goal of 
1% jet energy uncertainty would yield a systematic error of about 0.5% on the vector boson 
+ jets cross sections. 

For a Higgs boson mass Mh ^ 2Mw, the most promising Higgs discovery channel at the 
Tevatron and the LHC is its decay into vector boson pairs. In the decay into TF-bosons, 
no clear mass peak is observable, since the neutrino from the VF-decay leaves the detector 
unobserved. To establish a Higgs boson discovery in this channel, a precise understanding 
of the Standard Model background processes yielding vector boson pairs is mandatory. At 
present, vector boson pair production is described theoretically at NLO. At this order, vector 
boson pairs are produced from quark-antiquark annihilation and quark-gluon scattering. 
Gluon-gluon fusion into vector boson pairs contributes only at next-to-next-to-leading order 
(NNLO), but could yield a potentially large contribution at the LHC because of the large 
gluon luminosity. 

The NNLO gluon-gluon fusion contribution to vector boson pairs production was recently 
computed in [3] , including all vector boson decay information. Although gluon-gluon fusion 
contributes only four per cent to the total vector boson pair production cross section, its 
importance is substantially enhanced to above ten per cent by particle reconstruction cuts, 
and further to thirty per cent by Higgs boson search cuts. This observation clearly highlights 
the need for a full NNLO calculation of vector boson pair production. First steps in this 
direction have recently been completed with the calculation of the two-loop and one-loop 
squared corrections [4] to the quark-antiquark annihilation matrix elements in the high 
energy limit. 



3 Underlying Event 



All Three Regions Charged PT Sum Density: dPT/dqd^ | 



r P T > 0.5 GeV/cimd In,! <1 



8 



CDF Run 2 Preliminarx 



:orrected Toward Data 
Away Particle Level MC (PYTHIA Tune AW) 
Corrected Away Data 



ZRegion: 70 < M < 110 GeV/c 



t f 



! t 



Transverse Momentum of Lepton Pair (GeV/c) 



The term "underlying event" summarises all 
event activity which is observed besides the 
hard interaction, for example: multiple in- 
teractions, event pile-up and remnant inter- 
actions. Since this underlying activity re- 
sults typically only in low-momentum par- 
ticles, it is not accessible to a fully perturba- 
tive description (although at least multiple 
interactions may be calculable in the high 
energy limit [5]). At present, the description 
of the underlying event in Monte Carlo gen- 
erator programs is based on models, such as 
the eikonal model of multiple scattering [6] 

which is used in HERWIG+-I- [7]. These models require tuning to experimental data, and 
their extrapolation from Tevatron to LHC energies is highly uncertain. 

In order to study and fit the underlying event in data, regions can be defined where the 
underlying event contribution is enhanced. Figure [2] shows an example of dividing Tevatron 
Drell-Yan events into three regions with respect to the <j> direction of the lepton pair [8]: 
the "toward" direction is defined as a cone within A(f> < 60° to the Z boson, the "away" 
direction is defined as a cone within < 60° in the opposite direction, and the "transverse" 



Figure 2: Underlying event study in vector bo- 
son production at CDF. Figure taken from [5]. 
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direction, which is very sensitive to the underlying event, is defined as the remaining area. 
The three regions are compared to PYTHIA, tuned to a similar analysis using high px jet 
events (PYTHIA tune AW), and show good agreement. 

The impact of the underlying event will be much larger at the LHC, where one expects 
roughly 35 minimum bias events per bunch crossing at the design luminosity of 10 34 cm _2 s _1 . 
The CMS experiment presented plans [9] for measuring the minimum-bias contribution using 
a dedicated forward hadronic trigger and also presented studies of the remnant interactions 
besides the hard scattering, using the same analysis strategy as at the Tcvatron. These 
studies are part of an overall strategy to quickly tune and understand QCD reactions at the 
new energy frontier. 

4 Monte Carlo Tools 

Event generator programs are based on leading order matrix elements, describing a primary 
(low multiplicity) hard scattering process, which is then used as input to a parton-showcr 
to generate higher multiplicities. The resulting multi-parton final state is then transformed 
into a multi-hadron final state using a hadronisation model. Based on Monte Carlo algo- 
rithms, these programs provide samples of unweighted events, which can then be further 
processed using detector simulation and event reconstruction software. They are employed 
very widely in all aspects of experimental studies at particle colliders. The programs which 
are currently used most extensively were initially developed about twenty years ago and have 
undergone continuous upgrades. Nevertheless, these programs now start to display serious 
shortcomings, since a variety of new theoretical developments (for example, the matching 
of leading order multi-particle matrix elements onto the parton shower, or higher order cor- 
rections or improved shower prescriptions) cannot be incorporated into them. To overcome 
these limitations, several completely new Monte Carlo codes are currently being developed. 

The HERWIG++ project [10] is a new Monte Carlo generator program, aiming to incor- 
porate the ideas of the well-established HERWIG generator, which was most widely used in 
the previous generation of collider experiments. Its new fully functional release became avail- 
able recently, and it includes, among other improvements, standard interfaces to specialised 
matrix clement generators, simulation of a variety of beyond-Standard-Model reactions, a 
consistent treatment of radiation off heavy particles, and the simulation of underlying event 
dynamics 111] using an eikonal model [6]. 

The SHERPA project [12] is a newly developed Monte Carlo generator program, aiming 
to incorporate many of the recent new theoretical developments: recent additions in this 
project include the merging of multi-parton matrix elements with parton showers |13j . a 
new shower model based on the dipole formalism [14] , preparations for automated NLO 
calculations [15j . and a new matrix-element generator based on improved Bcrcnds-Gicle 
recursion relations [16] , 

All generic multi-purpose Monte Carlo programs are at present restricted to leading order 
in perturbative QCD. The extension of Monte Carlo programs to include NLO corrections 
is a currently ongoing activity, and has been accomplished for a variety of specific processes 
already [T7] . For internal consistency, leading order Monte Carlo programs should use parton 
distribution functions at leading order. In fitting LO parton distributions to deep inelastic 
scattering and hadron collider data, one observes a poor fit quality, largely due to the high 
precision of the experimental data, which arc sensitive to higher order QCD effects in the 
different observables. To improve the quality of the LO description of data on the proton 
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structure, various modifications were suggested [18] : by easing the momentum sum rule and 
modifying the scale in the QCD evolution, it is possible to mimic some of the numerically 
dominant higher order effects. An alternative approach emerging from discussions in the 
working group would be to consider the parton distribution functions to be an integral part 
of each Monte Carlo event generator, and to include them as parameters to be tuned. 

5 QCD in the High Energy Limit 

In the conventional fixed-order approach to perturbative QCD, scattering cross sections 
are computed as expansion in the strong coupling constant, and the structure of incoming 
hadrons is described by parton distributions obtained within collincar factorisation, evolving 
according to the DGLAP evolution equations. This fixed-order approach provides a very 
successful description of a broad range of observables, it becomes however inappropriate if 
higher-order terms in the coupling constant expansion are enhanced by large logarithmic 
corrections, which can spoil the convergence of the perturbative series. In these cases, 
an all-order resummation of the large logarithmic corrections is required to obtain reliable 
predictions. In the high energy limit of QCD scattering processes, which corresponds to low 
x in deep inelastic scattering, terms of the form cc™ ln m x can become potentially large at all 
orders, thereby invalidating the fixed-order expansion. In this limit, collinear factorisation, 
which assumes transverse-momentum ordering of initial state radiation, becomes equally 
inappropriate and should be replaced by transverse-momentum (kr) factorisation, yielding 
unintegrated parton distributions. The resummation of large logarithms in the high energy 
limit is accomplished by the BFKL evolution equation. 

Since the inclusive proton structure 
function at HERA energies is not suffi- 
ciently sensitive to differences between the 
DGLAP and BFKL approaches, experimen- 
tal studies of BFKL resummation effects fo- 
cus largely on specific hadronic final states 
at small x. Especially the forward jet cross 
section has turned out to be very discrimi- 
native between different approaches. 

Data on this observable, single, dou- 
ble and triple differential cross sections are 
in general only poorly described by NLO 
QCD. Among standard parton shower event 
generator programs, only ARIADNE agrees 
with observations [T9l [20] , after having been 
tuned to other HERA data. ARIADNE is 
based on the colour dipole model and ex- 
hibits BFKL- like parton showers unordered Fi g ure 3: Di J et cross sections normalised to 
in k T . Differential information may also be the visible cross section between 0° and 170° 
gained from azimuthal correlations of dijet at HERA compared to NLO 3-jet and NLO 
production at low x [S], which is not fully 2 ~i ct calculations. Figure taken from [21]. 
described by NLO QCD as illustrated in 

Figure [H This observable is potentially sensitive to the unintegrated gluon distribution [22] . 
and may enter into a global determination of unintegrated parton distribution functions. 
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The BFKL equation at leading logarithmic (LL) approximation correctly describes the 
overall features of forward jet production. A fully reliable description can however only 
be attempted by including subleading logarithmic corrections (NLL), which are at present 
not yet fully available for jet production in deep inelastic scattering. An approximate NLL 
study [23j does show clear improvements upon BFKL at LL, such as a stabilisation of 
the scale and scheme dependence and agreement with experimental data from HI over an 
extended kincmatical range. 

To compute cross sections within fey-factorisation, one needs to derive scattering ampli- 
tudes with off-shell initial state partons. To obtain those in a gauge-invariant form, it is 
most convenient to couple the off-shell partons to an external current [24] . Recent calcula- 
tions in this framework have focused in particular on final-state photon production |251 I26j 
in proton-proton and photon-proton collisions. Since these calculations do not provide the 
full final state information on all partons, issues like photon isolation and infrared cut-offs 
are still controversial. 



The BFKL evolution equation at leading 
logarithmic accuracy can be reformulated as 
a parton shower by interpreting its evolu- 
tion kernel as splitting probability [27]. In 
this formulation, constraints from energy- 
momentum conservation are implemented 
in a straightforward manner. The unin- 
tegrated parton distribution functions re- 
quired as initial conditions are inferred from 
ordinary DGLAP parton distributions by 
undoing the last branching in the DGLAP 
evolution. With this parton shower for- 
mulation, leading logarithmic BFKL pre- 
dictions can be made for a great variety 
of observables. For example, inclusive jet 
production at the Tevatron is correctly de- 
scribed over a wide kinematical range, Fig- 
ured] It can be clearly seen that in certain 
kinematical regions different jet multiplici- 
ties can be of comparable numerical magni- 
tude in the BFKL approach. 
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Figure 4: Inclusive jet cross sections at Teva- 
tron compared to LL BFKL parton shower 
predictions. Figure taken from [27| . 




The BFKL formalism can also be ap- 
plied directly to approximate the matrix elements for hard scattering needed for the predic- 
tion of multi-jet events. When supplemented by simple constraints on the analytic behaviour 
of the amplitudes [28] , a good approximation of the corresponding full leading order matrix 
elements is acheived, while allowing for an all-order resummation of the hard perturbative 
corrections. 

A complementary approach to the high-energy limit of QCD is the eikonal approximation, 
which allows to compute soft-gluon corrections to all orders in the coupling constant. This 
approximation is applied for example in pair production [29] of heavy quarks. 
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6 Hadron Spectroscopy 
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Hadron spectroscopy has been playing an 
important role in understanding strong in- 
teractions. Recently, ZEUS has stud- 
ied KgKg final states in ep collisions 
at HERA 30]. Three enhancements in 
invariant-mass distribution were observed, 
as shown in Figure El They correspond to 
/a(1270)/ag(1320), /a(1525) and /„(1710) 
mesons. The interference pattern, predicted 
by SU(3) symmetry, was taken into account 
fitting the invariant-mass distribution. The 
/o(1710) state, which has a mass consistent 
with a J PC = ++ glueball candidate, is ob- 
served with 5 standard-deviation statistical 
significance. However, if this state is the 
same as that seen in 77 — > Kg Kg [31] it is 
unlikely to be a pure glueball state. 

A broad spectrum of results on hadronic 
interactions at low energies, including nuclear form factors and resonances 
transparency [33] is currently being obtained at TJNAF. 
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7 Heavy-Ion Collisions 

At the Relativistic Heavy-Ion Collider, jets 
have been used to study the high energy 
density matter created in nuclear collisions. 
Defining regions in a similar way as in un- 
derlying event studies in pp collisions, one 
can study the jet evolution in the quark 
matter. Figure [34] compares the jet evolu- 
tion in pp and Au-Au collisions for increas- 
ing jet energy (here a jet is defined as a 
high energy hadron). In this figure, solid 
histograms (shaded bands) indicate ellip- 
tic flow model uncertainties. Arrows in (c) 
show the "head" region (HR) the "shoul- 
der" region (SR) and the "near-side" re- 
gion (NR) . The pp data shows a typical two 
peak structure due to back to back dijets. 
The Au-Au spectra reveal jet quenching in 
the medium as the pt increases, and in ad- 
dition, they show a more complicated jet 
structure evolution with px, with prominent 
tions - a feature expected from a medium-in 
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Figure 6: Per-trigger yield versus A<p for 
various trigger and partner pt arranged by 
increasing pair proxy energy. Figure taken 
from [M]. 

peaks in the "shoulder" regions at fixed posi- 
luced Mach shock. 
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8 Hadron Fragmentation 

Charged particle production in DIS at 
HERA was studied in terms of multiplicity 
and scaled momentum distributions [35(136]. 
When an appropriate energy scale is used, 
ep data can be consistently compared with 
data from e + e~ , fip and vp scattering over a 
wide energy range. In all cases, similar be- 
haviour is observed, supporting quark frag- 
mentation universality, see Figure [7] 

Models, implemented in LO matrix ele- 
ment Monte Carlo programs, describe data 
reasonably well but NLO QCD calculations, 
using three different fragmentation func- 
tions, fail to describe the scaling violations 
seen in the data [36] , see Figure [8] Figure 7 : Mean charged multiplicity. Figure 

The production of K S ,A and A in DIS taken from [35] . 
at HERA has been investigated by HI [37] , 

The model predictions, based on leading order Monte Carlo programs in general are able to 
describe the overall features of the measurements, however they fail in details, in particular 
in describing A and A production in the current region of the Breit frame. It was also found 
that a constant strangeness fraction in hadron fragmentation fails to fit all the data. 

Data from pp, e + e~ and heavy ion scat- 
tering has been used to test the limiting 
fragmentation hypothesis [38] . and a pre- 
diction has been made for ep DIS |39| . 

Bosc-Einstcin correlations of hadron 
pairs at HERA have been studied by HER- 
MES [40] for nuclear targets ranging from 
hydrogen to xenon. It was found that the 
parameters describing the correlations nei- 
ther depend on the nuclear target nor on 
the hadronic invariant-mass. 

Momentum distributions of identified 
hadrons can be used to extract hadron frag- 
mentation functions, which describe the dif- 
ferential probabilities for specific parton-to- Figure 8: Normalised distributions of the 
hadron transitions. These fragmentation scaled momentum. Figure taken from [36] . 
functions obey DGLAP evolution equations 

with timelikc splitting functions. At present, these timelike splitting functions are known 
to NLO, and calculations to NNLO are in progress [5T]. Using extensive data sets taken in 
e + e~, pp and pp collisions, and including hadron mass corrections, a new global extraction 
of parton fragmentation functions to n^, K^, p/p, Kg and A/A was performed recently [42] . 
Compared to earlier studies, especially the determination of non-singlet quark fragmenta- 
tion functions is improved by charge asymmetry pp data from RHIC. Unfortunately, data 
on momentum distributions of identified light hadrons from HERA have up to now only 
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been released by the fixed-target HERMES experiment, while HI and ZEUS measurements 
of charged and neutral light hadron production released up to now were not suitable for the 
extraction of fragmentation functions. 
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9 Isolated Photons 

Direct photons coming from the partonic 
hard scattering are a powerful probe of of 
the dynamics of hard QCD interactions. At 
Tevatron energies, the direct photon + jet 
cross section measurement is sensitive to 
the gluon density and thus can be used to 
complement the HERA data quark on the 
gluon at low x by constraining the gluon 
contribution at large x. In addition, pho- 
ton final states are predicted in new physics 
models such as SUSY, Extra Dimensions, 
Technicolour, etc. The DZero experiment 
has measured the isolated (hence en- 
hanced direct) 7 + jet inclusive cross sec- 
tion in four regions of jet and photon rapid- 
ity and compared them with theory using 
JETPHOX [44] and CTEQ6.5M PDFs as 
shown in Figure The two dotted lines 

show the effect of theoretical scale variations by a factor of two; the shaded region indicates 
the CTEQ6.5M PDF uncertainty, and the dashed and dashed-dotted lines show ratios of the 
JETPHOX predictions with MRST 2004, Alekhin and ZEUS 2005 PDF sets to CTEQ6.5M. 
The systematic uncertainties have large bin-to-bin correlations in p7p. An additional 7.8% 
normalisation uncertainty on the data points is not shown. 
While there is overall agreement with 

Inclusive prompt photons in photoproduction at H1 
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Figure 9: Ratio of measured and predicted 
cross section for pp — > j+jet+X as a function 
of p\ . Figure taken from [43] . 
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the theoretical prediction, and the theory 
successfully models the pr distributions for 
the isolated 7+jet measurements at HERA, 
the theory does not describe well the pTp 
dependence of the cross sections as mea- 
sured 'mpp at the Tevatron. Neither reason- 
able variations in fragmentation functions 
nor contributions from threshold resumma- 
tion are able to improve the data description 
from theory. Clearly this is an important 
measurement to be understood in prepara- 
tion for data taking at LHC. 

First results on prompt photons in pho- 
toproduction by the HI collaboration mak- 
ing use of the full HERA-II statistics were 

presented [55]. They cover the phase space 5 < < 15 GeV and —1.0 < rp < 2.4 in the 
laboratory frame. The isolation criterion for the prompt photon is defined by the rcquirc- 
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Figure 10: Inclusive prompt photons in photo- 
production at HERA. Figure taken from [45] ■ 
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merit that E^/E^T 3et > 0.9. The data are above the two QCD predictions, particularly 
at low Ej, as can be seen in Figure [TO] The FGH prediction [IS] is based on collinear fac- 
torisation at NLO and the ZL prediction on a fc^-factorisation approach [47] . They equally 
investigated prompt photon plus jet data are found to be reasonably well described by both 
predictions, except at the highest x° bs , where direct photoproduction is enhanced. 



10 Precision Physics with Jets 

Inclusive jet cross sections at pp colliders are a sensitive probe of parton distribution func- 
tions and of new physics. The DZcro experiment has made the most precise measurement 
of the inclusive jet cross section to date [48]. Figure [TT1 shows the ratio of the measured 
to predicted inclusive jet cross sections as a function of jet pt in six rapidity bins. The 
theoretical prediction comes from a NLO QCD calculation using the CTEQ6.5 PDFs. The 
MRST2004 PDF predictions are also plotted. As is demonstrated in this figure, the experi- 
mental uncertainties are now smaller than the scale and PDF uncertainties. The dominant 
experimental uncertainty comes from the jet energy scale, and after seven years of work, the 
DZcro collaboration has reduced this uncertainty to about 1.2% for high px central jets. 

Studies are ongoing at the LHC in preparation for data taking in late 2008. New physics 
can be seen as one goes to higher px in the inclusive jet spectrum. The current limit on the 
contact interaction scale at the Tevatron is 2.7 TeV. The CMS experiment at the LHC will 
be able [49] to improve this limit within the first 10 pb _1 of collected data. Assuming a jet 
energy scale uncertainty of 10% the predicted sensitivity to interaction scale is A ps 3 TeV, 
7 TeV and 10 TeV for 10 pb" 1 , 100 pb" 1 and 1 ftr 1 respectively. 
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Figure 11: Inclusive jet cross section as measured by DZero divided by theory as a function 
of pt for 6 \y\ bins. Measurement uncertainties are shown by the shaded bands and PDF 
uncertainties by the dashed lines. The theoretical prediction using the MRST2004 PDFs is 
shown by the dotted line. The hashed area at the bottom of the plots shows the uncertainty 
when varying the renormalisation and factorisation scales between pt/2 and 2px- Figure 
taken from [48] . 
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Recent measurements by the ZEUS col- 
laboration of 2,3 and 4 jets in photoproduc- 
tion were presented [50] . Compared to ear- 
lier 2-jet results, the jet transverse energy 
requirements (E 3 T etl{2) > 20(15) GcV), the 
range in pscudorapidity (—1 < rf ei < 3) 
and the luminosity have been significantly 
increased in order to obtain a dataset suit- 
able for testing and constraining both the 
photon and the proton PDFs. The sensitiv- 
ity to different PDFs was shown by compar- 
ing the data to different available parametri- 
sations. For the 3 and 4-jets analysis the se- 
lection of the jet phase space was different: 
> 6 GeV and \n> et \ < 2.4. In addi- 



E 3 T et 



tion, the data were divided into a low and 
high mass sample with 25 < M n j < 50 GeV 
and M n j > 50 GcV in order to look for ef- 
fects of multi-parton interactions or under- 
lying events and to allow for tests of multi- 
jet description in event generators. In gen- 
eral both models HERWIG and PYTHIA 
largely underestimate the rate of multi-jet 
events. At high M n j however, they provide 
a reasonable description of the shape of the 
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Figure 12: Photoproduction cross section of 
3-jets as a function of the 3-jet invariant-mass 
at HERA. Figure taken from [5U] , 



data, while low M„ 



and low x° bs 



are only described once multi-parton interactions are 



introduced. A leading order QCD calculation for 3-jcts, corrected for hadronisation effects 
and effects of multi-parton interactions, is compared to data in Figure [12l illustrating the 
need for a NLO calculation. 



For jets with sufficiently high E^ one 
may expect the contribution of fragmenta- 
tion to the jet substructure to become small 
compared to the contribution from parton 
radiation. Preliminary measurements by 
the ZEUS collaboration [51] of the mean 
integrated jet shape (^(r)) for two types 
of jets with 14 < E° T et < 17 GeV arc in 
good agreement with corresponding NLO 
calculations as illustrated in Figure [T31 The 
data are also well reproduced by the colour 
dipole model as implemented in ARIADNE. 
The "one-jet" distribution is mainly due to 
quark-initiated jets, selected by demanding 
events with 1-jet only. The "two-jet" distri- 
bution is enriched with gluon-initiatcd jets 
by choosing the jet with the lowest E^* in < 
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Figure 13: Integrated jet shapes for quark and 
gluon enriched jets in DIS at HERA. Figure 
taken from [5"T] . 

ents which have only two jets, which are close 



DIS 2008 



H1 Preliminary 



5 0.4 




HERAIPLB 653 (2007)134 

(phase space con'a} 

HERA Ml preliminary 
NLO®hadr 



to each other in rj — <j> space. One observes that the gluon enriched jets are broader as 
expected. 

Recent polarised and unpolarised inclu- 
sive jet cross sections were measured by 
the ZEUS collaboration in charged current 
(CC) e ± p events using HERA-II data [52] . 
They are found to be in good agreement 
with Standard Model predictions. These 
jet cross sections and the fully inclusive CC 
cross sections can provide constraints on the 
u and d-PDFs at high x. 

Preliminary results on normalised inclu- 
sive, 2-jet and 3-jet cross sections by the 
HI collaboration were presented |53| . using 
HERA-I and II data. In photon virtuality 
they cover the range 5 < Q 2 < 15000 GeV 2 
and in jet transverse energies 7(5) < E^* < 
50(80) GeV. The normalised inclusive jet 

cross section for the high <3 2 data, i.e. Q 2 > 150 GeV 2 , is shown in Figure [Ml together 
with published HERA-I data and NLO QCD calculations |33], illustrating good agreement 
and experimental errors, which are clearly smaller than the uncertainty of the NLO calcu- 
lation. These high-Q 2 data have reached an experimental uncertainty of about 3 %, which 
is mainly due to the uncertainty of the jet energy scale of about 1.5 %. 
The value of the strong coupling 
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Figure 14: Normalised inclusive jet cross sec- 
tion and ratio of data over theory at HERA. 
Figure taken from [5"3"j . 
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was extracted from a fit of the NLO calcula- 
tion to the normalised inclusive, 2 and 3-jet 
cross sections, resulting in an impressively 
small experimental error of 0.7 %. The total 
uncertainty is dominated by the uncertainty 
of the NLO theory, indicating the need for 
a NNLO calculation. The running of the 
strong coupling is verified over two orders 
of magnitude in Q as demonstrated in Fig- 
ure nu 

Using the recently computed NNLO cor- 
rections to event shape variables [55] , a new 
extraction of a s from data on the standard 

set of six event shape variables, measured [55] by the ALEPH collaboration at LEP1 and 
LEP2 was performed. One observes a clear improvement in the fit quality when going to 
NNLO accuracy [57]- Compared to NLO the value of a s is lowered by about 10%, but 
still higher than for NLO matched with the resummed next-to-leading logarithmic approx- 
imation (NLO+NLLA) (56]. As illustrated in Figure [TBI the scatter among the a s -values 
extracted from different shape variables is lowered considerably, and the theoretical uncer- 
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Figure 15: Running of a s from jet cross sec- 
tions at HERA. Figure taken from [53] . 
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tainty is decreased by a factor 2 (1.3) compared to NLO (NLO+NLLA). The combination 
of all shape variables at all energies yields 

a s (M z ) = 0.1240 ± 0.0008 (stat) ± 0.0010 (exp) ± 0.0011 (had) ± 0.0029 (theo). 
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The fixed-order QCD description of 
event shape distributions is reliable only if 
the event shape variable is sufficiently far 
away from its two-jet limit. In the approach 
to this limit, event shapes display large in- 
frared logarithms at all orders in pertur- 
bation theory, such that the expansion in 
the strong coupling constant fails to con- 
verge. Rcsummation of these logarithms 
yields a description appropriate to the two- 
jet limit. To explain event shape distri- 
butions over their full kincmatical range, 
both descriptions need to be matched onto 
each other. The matching of NLLA with 
NNLO was performed recently 58J. The 
most striking observation is that the dif- 
ference between NLLA+NNLO and NNLO 
is largely restricted to the two-jet region, 
while NLLA+NLO differ in normalisation 

throughout the full kincmatical range. This behaviour may serve as a first indication for the 
numerical smallncss of corrections beyond NNLO in the three-jet region. Fits of a s based 
on NNLO+NLLA are currently in progress. 

A summary of the most precise a s (Mz) 
extractions from HERA data discussed at 
DIS 2008 is shown in Figure [TTJ It also in- 
cludes the "HERA comb. 2007 incl. jets" 
result |59j of a first simultaneous fit to pub- 
lished HERAT inclusive jet cross section 
data from HI and ZEUS. 



Figure 16: Strong coupling constant extracted 
from ALEPH data on event shapes. Figure 
taken from [571. 
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11 New Jet Algorithms 

The reconstruction of hadronic jets at col- 
liders proceeds through the application of a 
jet algorithm, which clusters the detector- 
level information on the hadrons (hadronic 
tracks or calorimeter energy deposits) into 
a limited number of composite objects: the 
jets. One distinguishes two types of jet al- 
gorithms: cone-based and cluster-based al- 
gorithms. The cone-based algorithms aim 

to maximise the hadronic energy inside a cone of fixed size, while the cluster-based algo- 
rithms perform a sequential recombination of parton pairs. For both classes of algorithms 



Figure 17: Summary of extractions of a s (Mz) 
at DIS 2008. 
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many different realisations have been proposed in the past. Depending on the observable 
under consideration, one or the other type of algorithm may be more appropriate: while 
cone-type algorithms perform generally better in the reconstruction of resonances, i.e. in 
searches for massive particles decaying into jets, cluster-type algorithms are more appro- 
priate for precision studies. To compare experimental jet measurements with perturbative 
QCD, any algorithm of either type must fulfil infrared safety criteria. 

Unfortunately, most of the cone-based algorithms used up to now at the Tevatron did not 
fulfil those criteria and display infrared sensitivity above a certain final-state jet multiplicity. 
The principal cause of this infrared sensitivity is the use of cone seeds to speed up the jet 
reconstruction. The seedless infrared-safe cone algorithm (SISCone) [60] overcomes these 
problems and provides a cone-type algorithm for hadron collider physics. 

The sequential recombination used in clustering-based algorithms is very time-consuming, 
since the distance measure for each pair of (pseudo-)particles has to be evaluated in each 
iteration. The practical applicability of these algorithms was therefore severely restricted, es- 
pecially at hadron colliders, where typical events contain a very high multiplicity of hadronic 
objects. Using techniques from computational geometry, the clustering can be performed 
in a much more efficient way |61j . such that cone-based and clustering algorithms display 
similar performance |62j . 

By generalising the fcr-clustering algorithm to allow arbitrary powers in the transverse 
momentum weight, one can define a one-parameter family of clustering algorithms. The 
weighting with the inverse power of the transverse momentum defines the anti-fey algo- 
rithm 63J. A new powerful tool to analyse the features of jets is their catchment area [64] . 
which is obtained by including zero-energy ghost particles in the clustering. By inspecting 
the jet areas, one observes that the anti-fey algorithm yields perfect cones, which makes this 
algorithm an ideal replacement for iterative cone- type algorithms. 

The jet area may be turned into a pow- 
erful tool to disentangle the hard scattering 
process from any underlying event activity. 
On an event-by-event basis, this underly- 
ing activity, which appears to be uniform 
in rapidity, could be measured in areas out- 
side any jet, and then subtracted from the 
event prior to reconstruction. To illustrate 
the improvements which can be obtained us- 
ing this procedure, Figure [18] displays the 
effect of event pileup on the reconstruction 
of the new heavy Z' gauge boson decaying 
into hadronic jets. While the mass peak is 
substantially shifted by the effect of pileup, 
it is observed at its true mass after jet area 
based subtraction. 



0.015 



S 0-01 



0.005 



k t , R=0.7 

LHC, high lumi I 
Z' at 2 TeV 



no pileup 
no pileup, sub 
pileup 
pileup, sub 




2100 



Figure 18: Effect of jet area based subtrac- 
tion on reconstruction of heavy resonance at 
M = 2 TeV decaying into jets. Figure taken 
from l64l. 



12 Conclusions 



The broad spectrum of new results on "Hadronic final states and QCD" presented in the 
working group, as well as the very lively discussions, illustrate the importance of this subject 
area in the current transition period between HERA and the LHC. Many recent experimental 
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results are of direct relevance to preparations for LHC physics, and theoretical physics is 
responding to a number of challenges in complexity and precision posed by the physics 
programme of present and future colliders. 
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